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Abstract 

We revisit the issue of superconductivity at the quantum-critical point between a 2D paramag- 
net and a spin-density- wave (SDW) metal with ordering momentum (vr,7r). This problem is highly 
non-trivial because the system at criticality displays a non-Fermi liquid behavior and because the 
effective coupling constant A for the pairing is generally of order one, even when the actual interac- 
^ ■ tion is smaller than fermionic bandwidth. Previous study [M. A. Metlitski, S. Sachdev, Phys.Rev.B 

o; „ 

^ ■ 82, 075128 (2010)] has found that the leading renormalization of the pairing vertex contains log , 

like in color superconductivity. We analyze the full gap equation and argue that summing up log^ 
term does not lead to a pairing instability. Yet, superconductivity has no threshold, even if A is 
set to be small: the subleading log terms give rise to BCS-like oc e-^A. We argue that the 
analogy with BCS is not accidental as superconductivity at a QCP is a Fermi liquid phenomenon 
- it comes from fermions which retain Fermi liquid behavior at criticality. We computed Tc for the 
actual A and found a good agreement with the numerical results. 
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Introduction. Superconductivity at the onset of density-wave order in a metal is an 
issue of high current interest, with examples ranging from cuprates jl], to Fe-pnictides ^ 
and other correlated materials 3|-|5| It is widely believed that the pairing in these systems 
is caused by repulsive electron-electron interaction, enhanced in a particular spin or charge 
channel, which becomes critical at the quantum-critical point (QCP). The pairing problem 
at QCP is highly non-trivial in D < 3, as scattering by a critical collective mode destroys 
Fermi liquid (FL) behavior above (Ref. js ol). This is particularly relevant for systems 
near uniform density- wave instability (e.g., a ferromagnetic or a nematic one). In this case, 
FL behavior is lost on the whole Fermi-surface (FS), and superconductivity can be viewed as 
a pairing of incoherent fermions which exchange quanta of gapless collective bosons j7-12|. 



The pairing of incoherent fermions is qualitatively different from BCS/Eliashberg pairing of 
coherent fermions in a FL because in the incoherent case the pairing in D < 3 occurs only 
if the interaction exceeds a certain threshold jsl lislll^ . For D = 3 there is no threshold, but 
at small coupling constant A, log A/Tc = l/\/X rather than 1/A (Ref. 16|), in close an alog y 
to Tc in color superconductivity (CSC) of quarks mediated by the exchange of gluons [17 1 

The non-FL behavior at criticality is less pronounced for systems near density-wave order 
at a finite momentum, because only fermions near particular points along the FS (hot 
spots) lose FL behavior at criticality. Still, fermions from hot regions mostly contribute to 
the pairing, and early studies of superconductivity at the onset of (tt, tt) spin-density-wave 
(SDW) order placed the pairing problem into the same universality class as for QCP 



with g = 0. The 2D problem has been recently re- analyzed [iQj] by Metlitski and Sachdev 
(MS). They argued that it is important to include into the consideration the momentum 
dependence of the self-energy along the FS, neglected in earlier studies. Using the full form 
of S(ci;m,k) for k on the FS, they found that the one-loop renormalization of the pairing 
vertex is larger than previously thought - it is log^ instead of log, and that the enhancement 
comes from fermions somewhat away from hot spots, for which E(co'm,k) has a FL form 
at the smallest frequencies. The log^ behavior in the perturbation theory holds for CSC, 
and MS result raises the question whether the pairing problem at a 2D SDW QCP is in 
the same universality class as CSC. The related issues raised by MS work are: (i) is the 
problem analogous to the pairing at a 2D SDW QCP a FL phenomenon, or non-FL physics 
is essential, (ii) what sets the scale of T^, and (iii) is Tc non-zero only if the coupling A exceed 
a finite threshold, as it happens if one approximates T,{uJm, k) by ^(wm) at a hot spot, or Tc 
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is non-zero even at smallest A, like in CSC? 

In this letter, we address these issues. We first show that the analogy with CSC does not 
extend beyond one-loop order, and in our case the summation of log^ terms in the Cooper 
channel does not give rise to a pairing instability. However, that subleading log terms do 
give rise to a pairing instability, and at weak coupling yield logA/T^ = 1/A, like in BCS 
theory. We show that the analogy with BCS formula is not accidental because the pairing 
predominantly comes from fermions for which fermionic self-energy has a FL form. In this 
respect, non-FL physics at a QCP is actually irrelevant to the pairing. We then analyze 
the physical case A = 0(1) in the large N approximation, assuming that terms with n > 6 



loops, which break expansion 18|jl9|j2l|, do not change the analysis qualitatively. We 



argue that the scale of is set by the frequency loq at which T,{ijjm) at a hot spot becomes 
equal to Um- The theory is under control when uq < Vp/a {a is inter-atomic spacing), which 
in practice implies that the interaction must be smaller than the fermionic bandwidth. We 
show that the prefactor c in Tc = cojo depends on the geometry of the FS and is c = 0.04 
when the velocities at the two hot spots separated by (tt, tt) are nearly orthogonal. This 
number agrees with the slope of Tc obtained by solving the gap equation numerically along 



the full FS 



22|. 



The model. We follow earlier works |8lJ9lJl9| and consider a system of fermions 
with N flavors and a large FS, like in the cuprates. Fermions are interacting by exchanging 
quanta of collective excitations in the spin channel. The static propagator of collective spin 
excitations is created by fermions with energies of order bandwidth and is an input for 
the low-energy theory of superconductivity. We assume that collective spin excitations are 
peaked at Q = (vr, vr) and focus on the hot regions on the FS, for which k^? + Q is also near 



the FS. The Lagrangian of the model is given by 

-A 



19| 



+9 / 'ipi+g,a^al3'ipk,l3 " S_g. (l) 
Jk,q 

where stands for the integral over rf— dimensional k (up to some upper cutoff A) and 
the sum over fermionic and bosonic Matsubara frequencies, Gq {k) = Go{um, k) = l/i^ium — 
VF,k(k — k^)) is the bare fermion propagator, and Xo (?) = Xoi^m, q) = Xo/(ct^ + ^^^) is the 
static propagator of collective bosons, in which measures a distance to a QCP and q is 
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measured with respect to Q. We set = below. The fermion-boson couphng g and xo 
appear in theory only in combination g = g'^Xo and we will use g below. The Fermi velocities 
at hot spots separated by Q can be expressed as = {vx,Vy) and Yp^2 = {—Vx,Vy), where 
X axis is along Q. We will also use a = Vy/v^ and vp = {v^ + VyY^"^. The model of Eq. 
([T]) can be equivalently viewed as a four-patch model for fermions near hot spots at ±k^ 
and ±(kir + Q) (Ref. 19|j2l|). The hot spot model is obviously justified only when the 
interaction g is smaller than Ep. 

The fermion-boson coupling gives rise to fermionic and bosonic self-energies. In the 
normal state, bosonic self-energy accounts for Landau damping of spin excitations, while 
fermionic self-energy accounts for the mass renormalization and a finite lifetime of a fermion. 
At one-loop level, self-consistent normal-state analysis yields 



x(^m,q) 



Xo 



q2 + \n^\'y 
3^ 



(2) 



Anvf 



l\uJm\ + 



2fc||a 
1+^ 



2k\\a 



(3) 



where 7 = 2Ng/{nVxVy) and fcy is a deviation from a hot spot along the FS. Right at a hot 
spot S(wm,0) = Um\/oJo/\oJm\, wheie Uq = {9g/{lQTTN)){2vxVy/vp). In ([2]), (|3]), we neglected 
terms with n > 2 loops, which are small in and terms with n > Q which break 
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expansion 
when = 1 (Refs. 
numerically small. 



19|. Superconducting is only weakly affected by two-loop corrections, even 



231]), and we expect that corrections from higher loops will remain 



Pairing vertex We add to the action the anomalous term $o(fc)V'fc,a(^o'^)a/3'i^-fc,/3 
and use Eq. ([1]) to renormalize it into the full $(A;). At T^, the pairing susceptibility 
Xpp{k) = ^{k)/^o must diverge for all k. The bare $0 can be set constant within a patch, 
but has to change signs between patches separated by Q (the pairing symmetry at the onset 



of SDW order is c?— wave 
obtained by MS: 



24|). The one- loop renormalization of $(A;) at k 



[UJ 



$(cc;~r,0) = <l'o(l + Alog^A/T), A 



a 



T, 0)was 



(4) 



7r(l + a2)' 

where A is the smaller of uq and a^Ep/uQ. Notice that neither the coupling constant g 
not 1/A^ are present in (j4]), the only parameter is the ratio of the velocities a, which is 



a geometrical property of the FS. For a cuprate-like FS, a ~ 1, i.e., the pairing couphng 
constant A = 0(1). To understand the physics of the pairing at the QCP, it is, however, 
advantageous to keep A as a free parameter and analyze the pairing at both small and 
moderate A. 




FIG. 1: Diagrammatic representation for the pairing vertex. The shaded triangle is the full 
the unshaded vertex is the bare ^o, solid lines are full fermionic propagators, and the wavy line 
is the Landau-overdamped spin propagator. The pairing vertex contains io"^^, the vertices where 
wavy and solid lins meet contain a^s- 

Let's first see where log^ renormalization comes from. The one- loop diagram for $ con- 
tains two fermionic propagators G{k) and G{—k) and and one bosonic x(^) (FigH])- Large 

allows one to restrict xi^m, k) to momenta connecting points at the FS and integrate 
over momenta transverse to the FS in the fermionic propagators only. Because S does 
not depend on this momentum, the integration is straightforward, and yields, to logarith- 
mic accuracy / GGx oc f dk\\ frj.dQrn{x{^m,k\\)/\^m + ^{^m, k\\))\. At kj > 'jQm and 
|A;||| < kpg/vp, + scales as and oc Integrating 

over A;|j we obtain J^^^ | dk'^^/k'^^ oc log |fim|, and the remaining integral over frequency yields 
/ GGx oc Jrp{dQ,rn/\^m\) log|i^m| OC log^T. We see that the log^ T dependence originates 
from extra logarithm from fc— integration. This in turn is the consequence of Qrn/k\\ form of 
self-energy ^||) at /cy > 7^7^- As S oc is the property of a FL, the log^T renormal- 

ization comes from fermions which preserve a FL behavior at a QCP. We further see that the 
one-loop renormalization can be interpreted as coming from the process in which fermions 
are exchanging quanta of an effective local log Q interaction. The same process determines 
one-loop renormalization of $ in CSC. 

The log^ analysis can be extended beyond leading order. We assume that A is small 
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FIG. 2: Numerical solution of Eq. ^ at small A. (a) The transition temperature. When A 
decreases, InXloguo/Tc approaches 1, as in Eq. ([8]). (b) The eigenfunction ^{y), where y = 
/cy / {ttT^) . Solid and dashed lines are numerical and analytical solutions of Eq. [6l respectively. 
The two are very close, except for the largest y ~ ojq/T, when the cutoff becomes relevant. 

but Alog^T = 0(1), and sum up ladder series of A log^ T terms neglecting smaller terms at 
each order of loop expansion. Performing the calculations (see Supplementary material for 
details), we find that the the analogy with CSC does not extend beyond leading order: for 
CSC the summation of A log^ T terms yields $ = $o/ cos[(2Alog^ rj^/^] (Ref.pjj), and the 
system develops a pairing instability at | logT^ = 7r/2v^2A (Ref. 17|). In our case, pertur- 
bation series yield $ = <l>oe^'°s -^ ^^ 

pairing susceptibility increases with decreasing 

T, but never diverges. 

Because the summation of the leading logarithms does not lead to a finite one has to 
go beyond the leading logarithmical approximation and analyze the full equation for $(fc) 
at $0 = in order to understand whether or not Tc is finite at a QCP. This is what we do 
next. 

Full gap equation. Within our approximation, the full linearized equation for the 
anomalous vertex is obtained by summing up ladder diagrams and keeping the self-energy 
in the fermionic propagator. Integrating the r.h.s. of this equation over momenta transverse 
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to the FS, we obtain 

3^ f dk\^ ^{um',k\^) 



$(a;^,/c||) = — T > / — — -7 



fcy + k'^ — 2ek\\k'^^ + 'y\ujm — oj^i \ 

where e = (1 — q;^)/(1 + a^). The temperature at which the solution exists is T^. Let's 
assume first that A is small and search whether is non-zero. At small A, typical fcy are 
larger than typical '^oJmi and one can easily verify that the vertex (^{ujmik\\) has a stronger 
dependence on k\\ than on frequency. In this situation, one can approximate (^{ujmik\\) by 
explicitly sum up over frequency and reduce ([5]) to ID integral equation in momentum. 
Introducing T = vrT/cJo and x = k'^^/{'yujQT) and setting a = 1 for simplicity, we obtain from 

m 

= A / (6) 



Ji 

The term in denominator with V xT is a soft upper cutoff. 

The r.h.s. of ([6]) contains log^ contributions from the range x ^ y, but, as we just found, 
they do not lead to a pairing instability. We therefore focus on the contribution from x ~ 
Because the kernel is logarithmical, we search for $(a;) in the form = exp[— /(p(a;))], 
where p{x) = 27rAlogx. Substituting this into we find that the form is reproduced at 
1 ^ X ^ l/T, when soft cutoff can be omitted. The self-consistency condition yields (see 
Supplementary material) 



fiz) = arcsin z + Vl - z'^ - Ij ■ (7) 

At small A, the soft cutoff can be replaced by the boundary condition that df{z)/dz must be 
at a maximum at z = 27rA| logT| (i.e., that $(y) rapidly drops at x ~ I/?)- This condition 
selects T = Tc to be 

T, ~ uoC'/^'-'^ (8) 

To verify this reasoning, we solved Eq. (jG]) numerically. We found very good agreement with 
analytical results (see Figj2]). 

We see from ([8]) that is non-zero already at infinitesimally small A, like in BCS theory. 
The analogy is not accidental as the pairing predominantly comes from momenta away 
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from hot spots, for which a; ~ T, i.e., ~ ['~^ujqY^'^ ~ kpg/Ep. For g <^ Ep, when the 
calculations are under control, maximal k\\ is still much smaller than kp, i.e., the pairing is 
confined to a region near a hot spot. At the same time, typical 70; > 'yT^ are much smaller 
then jcoq, fermionic self-energy for A;|| ~ ('juo)^^'^ has the FL form {T.(ujm,k^\) oc uj„i/\kii\), 
and for x ~ ?/ in (|6]) the momentum integration does not give rise to an additional logarithm. 
The instability at Tc is then a conventional Cooper instability of a FL with a weak and non- 
singular attractive interaction. In other words, for small A, the pairing at a SDW QCP is 
entirely a FL phenomenon. 

A non-zero at small A is the consequence of the dependence of the self-energy on 



the momenta along the FS. Earlier works 



9| neglected this momentum dependence and 



approximated the self-energy by its non-FL form = uJm{^Q/\ujm\) ' at a hot spot. 

These studies found a different result: Tr becomes non-zero only if the coupling exceeds a 

nn 

certain threshold, like in the pairing problem at a QCP at g = (Refs. |14jj25l |. Specifically, 
because S = T,{um), the anomalous vertex <l> also depends only on frequency, and Eq. ([5]) 
reduces to ID integral equation in frequency rather than in momentum: 

iM = .'XTj: (9) 



where Z^^, = 1 + a/Iw^'I/wo- This equation has been solved for arbitrary A |8[, and the 
result is that Tc becomes non-zero only when A exceeds a critical value Ac = 0.035. Near 
critical coupling ~ cooe'^'^^^^^"^"^^^^), and for physical A = l/2n, = O.ITcjq- Note that 
the pairing in this case comes from fermions with non-FL self-energy, yet the overall scale 
for Tc is the same loq as in Eq. ([8]). 

Tc at moderate coupling. For a generic 2D FS, and Vy are comparable, and A 
from (jlj) is of order one. We solved Eq. Ofor Tc numerically for a = 1, when A = l/(27r), 
and found 

Tc ^ 0.04a;o. (10) 

Typical fcy ~ 'juq and typical 70; ~ 7Tc are now comparable, i.e., the pairing comes from 
fermions for which self-energy is in a grey area between a FL and a non-FL. Numerically 
Tc is still small compared to Uq, i.e., typical fcy are still larger than typical 'ju. To a good 
accuracy the self-energy of typical fermions can still be approximated by the FL form, i.e., 
even at A = 0(1) the pairing at a SDW QCP is still predominantly a FL phenomenon. 
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FIG. 3: Numerical solution of Eq. ([5]) for various a = Vy/vx (the Fermi velocities at hot spots 
separated by (7r,7r) are vf,i = {vx,Vy) and wf,2 = {—Vx-,Vy). We set ga/vp ~ 0.032. For a = 0(1) 
the momentum dependence of the fermionic self-energy dominates, pairing is a FL phenomenon, 
and Tc ~ 0.04a;o. At smaller a, when g > a^^'^vp/a (to the left of vertical dashed line), frequency 
dependence of the self-energy prevails, the pairing comes from fermions with non-FL self-energy, 
and Tc ~ O.ITwq, up to corrections of order {ga/vpY^'^ ■ In our case, these corrections reduce this 
number from 0.17 toto 0.15. 

Eq. ( [T0|) agrees with the numerical solution of the gap equation along the full FS, without 
restriction to hot spots [2^. In notations of Ref. |22|, Tc = vp/afiu), where dimensionless 
u = iuoa/^Svp)- Eq. ( ITO|) implies that f{u) = 0.03m at small u. This agrees well with the 



numerical solution. According to various numerics [26|, at larger u > 1/2, f{u) saturates at 
a value around 0.015 — 0.02 and at even larger u decreases as 1/m because of Mott physics. 
At strong nesting, when a -C I and 7 oc l/a is large, the pairing eventually become 
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determined by fermions whose self-energy has a non-FL form. This happens when {•yojoY^'^ 
becomes larger than maximal possible \k\\\ along a Fermi arc, which is of order 1/a, i.e., 
when vp/a ^ g > a^^'^vp/a. In this situation, the momentum dependence of T,{ujm,k\\) 
becomes irrelevant at Um ~ i^o; ^{^m,ki\) becomes predominantly frequency-dependent, 
and Tc recovers the value O.ITwq, which is for momentum-independent self-energy. We 
verified numerically that this is indeed the case. We plot Tc vs a in Fig. El At a = 0(1), 
Tc ~ 0.04(^0) and at small enough a approaches O.OITwq. We caution, however, that the 
limit a <^ 1 has to be taken with nesting may generate additional singularities at 

higher-loop orders. 

Conclusions. In this paper we analyzed the equation for superconducting at the 
onset of SDW order. We demonstrated that the leading perturbation correction to the bare 
pairing vertex contains log^ T, but the series of log^ T renormalizations do not give rise to 
the pairing instability. Yet, is finite, even when coupling is small, because of subleading, 
log T terms. We showed that the pairing predominantly comes from fermions which retain 
a FL form of the self-energy at a QCP, i.e., the pairing at a SDW QCP is actually a FL 
phenomenon. The overall scale of is set by the interaction (cuq ~ g), as long as the 
interaction is smaller than the Fermi energy. 

The large N approximation which we used to justify the gap equation requires caution 
because (i) two- loop 1/A^ corrections are logarithmically singular |9|j23|, and (ii) diagrams 



with n > Q loops are not small in 1/A^ and may give rise to additional logarithms |18lJl9lJ21 |. 
On the other hand, the fact that our Tc is of order uq, which is the upper limit of logarithmical 
behavior, likely implies that logarithmical renormalizations, not included into our analysis, 
do not change Tc qualitatively and do not affect the conclusion that the pairing at SDW 
QCP is a FL phenomenon. 

We acknowledge stimulating discussions with Y.B. Kim, S.S. Lee, M.A. Metlitski, S. 
Sachdev, T. Senthil, and A-M Tremblay. The research has been supported by DOE DE- 
FG02-ER46900. 
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I. SUPPLEMENTARY MATERIAL 



A. Perturbation Theory 

We add a fictitious anomalous term ^o4'k,a{io'^)ai34'-k,i3 to the original Lagrangian to 
generate a bare pairing vertex, renormalize it by particle-particle interaction, and obtain the 
pairing susceptibility as a ratio of the fully renormalized and bare pairing vertices. Within 
the approximations which we discuss in the main text, the diagrams for the renormalization 
of the pairing vertex form ladder series, shown in Fig. HJ however each line is the full 
fermionic Green's function, which includes the self-energy. 




FIG. 4: Perturbation expansion for the pairing vertex. 



The renormalized vertex $ depends on fermionic frequency u^, fermionic momentum 
along the Fermi surface, fcy, and on temperature. At high enough temperatures ^{uJm, k\\,T) 
weakly deviates from $o, but at Tc the ratio $/$o should diverge for all frequencies and 
momenta. For definiteness, we set k\\ to zero and Um to ttT, i.e., consider $(7rT, 0, T) = $(T). 
To simplify the formulas, we also set the ratio of Fermi velocities a = Vy/v^ to one. 

Consider one-loop renormalization of $ (the diagram with one wavy line in Fig. H]). 
After integration over momenta transverse to the Fermi surface (FS) we obtain, replacing 
the summation over Um by integration over \(jJm\ > t^T, 

dujmdk\\ 1 1 



$(T) = $0 1 + vr(a;o7) 



1/2 



47r2 \Urn + ^{0Jrn,k\\)\k'^\+'l\^i 



where, we remind, T,{ujm,k\\) = sign Um\'^{^^m, k\\)\, and [QIJIQI 

\T.{uJm,kii)\ = y/uJ^ (^1^1 + kf\/l - \h\/V^) ■ 
Rescaling = (|a;|7)^/^2; we obtain 
/ 

$(T) = $0 1 + A / ^= / dz 



1 



1 + 



\ 1/2 , 

j iVTT^+\z\)J 



(11) 



(12) 



(13) 
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where A = l/(27r). One can easily make sure that the the integral in the r.h.s. of (1T3|) 
contains \og^ uq/T, which comes from large z, i.e., from ^ ^ \uJm\l (Ref. 19|), and \oguo/T, 
which comes from z = 0(1). The coupling constant A is a fixed number, but in the spirit of 
perturbation theory we assume that A is a small parameter. Then, obviously, Xlog^ coq/T is 
more relevant than Alogwo/T, i.e., to logarithmic accuracy 19| 

$(T)<l>o(l + Alog'a;o/T) (14) 

We now check whether the series of Alog^Wo/T give rise to a pairing instability at 
\oguJo/T ~ l/VX. Because log^uJo/T comes from /cy ^ |i^m|7, we can expand the self- 
energy in 7|a;|//cy, i.e., approximate |S(a;m,, /cy)! in f|T2l) by 

|S(wm,A;||)| ~ (15) 
Using this self-energy, we obtain at two-loop order (/cy i = ki, /cy 2 = k2) 

$(T) = $o l + Alog^co/T + A^/ -1/ -^ + 0{\'^ 

(16) 

Evaluating the integrals we find 

$(r) = $0 (^1 + A log' tuo/T + log^ Uo/T + (17) 

To understand what are the prefactors from higher-order terms, we note that log^ term 
comes from the region of internal momenta k2 ^ ki, or, more specifically, the momentum 
ki in the cross-section which is farther from the vertex sets the lower cutoff of momentum 
integration over ^2 in the cross-section, which is closer to the vertex. This sheds light on 
the general pattern - we verified that at order M the largest, log^^^Wo/^ term comes from 
comes kM > ^m-i > ••• > ki, where kM is the momentum in the cross-section, which is the 
closest to the vertex. The same trend was earlier detected in the perturbation theory for tjbe 
pairing vertex in 2D systems for which fermionic self-energy depends only on frequency 
A simple exercise in combinatorics then yields, for the full series 

$(r) = (^1 + A log^ uo/T + ^A^ log^ uo/T + ^X' log' coo/T + ... 

M=0 
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We see that the calculation to log^ accuracy does not give rise to a pairing instability: the 
pairing susceptibility, $/$o increases with decreasing T, but does not diverge at any finite 
T. 



1. A comparison with color superconductivity 

The presence of \o^ujq/T in perturbation theory brings in the comparison with the 
problem of color superconductivity (CSC), where perturbative expansion also contains series 
of A log^ T terms. For that problem, however, previous works have demonstrated that the 
summation of log^ terms does lead to a pairing instability. It is therefore instructive to 
compare perturbation theory for CSC problem with our case and see where the two cases 
differ. 

In CSC problem 16lJl7l |. as well as in condensed-matter problems of the pairing medi- 
ated by gap less collective excitations in three spatial dimensions [lO|, fermionic self-energy 
depends only on frequency (and is actually irrelevant to the pairing problem), and the extra 
logarithm, in addition to a Cooper one, appears because the effective interaction, integrated 
over momenta along the FS, has logarithmic dependence on frequency. Because the integra- 
tion over momenta k\\ can be done independently in any cross-section in Fig. Hj the pairing 
can be analyzed within the effective local model of fermions with dynamical interaction 17 1 



x(wm) = Alog^^ (19) 

Substituting this form of the interaction (wavy line) into the one-loop diagram for $, we 
obtain 

$(T) = $0 (l + 2A £ ^ log ^) = (l + A log^ (20) 

This is the same result as in our case. At two-loop order, however, the prefactor for 
A^log'^A/T term is different from the one in our case. For CSC we have, at two- loop 
order 

$(T)=$o l + Alog2 — + AM -^\og— -^\og— 21 

V ttT J^t w k I At ^ \u - u \ J 

It is natural to divide the integration range over u" into two regimes, u" ^ u' and u" <^ u'. 
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It turns out that both regimes give contributions of order log A/T. We have 



2 :^2duj', A f^2duj'\ A 



J^T U' \UJ I J^T ^ \UJ -UJ \ 

,2f^ f^'^duj' A, A [""2(100' A, 2 A 

=AM 2 / ^ log ^ log — log - + / log ^ log' - 

X2 A 1 4 A 1 , 4 A 

=^A2log4A, (22) 
Combining one-loop and two-loop terms, we obtain 

$(T) = $0 (l + A log' A + i^g4 (23) 

which is different from the two-loop result in our case, Eq. f[T7|) . The difference comes about 
because for CSC there is no requirement that highest power of the logarithm comes from 
particular hierarchy of running frequencies - at two loop order, there is log^ contribution 
from the range where the highest frequency in the cross-section next to the vertex, and from 
the range when the highest frequency is in the cross-section farthest from the vertex. 

The series of A log' terms for CSC problem have been summed up in Ref. [igI, an the 
result is 

= ttT) = $0 f 1 + a log' A + i^g4 A ^ 

\ 111 b TVl 



$0 



cos[(2Alog'A/7rr)i/2] 
B. Linearized gap equation 

1. Interplay between characteristic momenta and frequency 



(24) 



,„ — , r dk',, (^(um',k'n) 1 
$(c.„,fc|,) = vr(c.o7)^/'T5^ 1 / "\ X . (25) 



We start with the Eq (5) in the main text. For a = 1 we have 

dk\^ ^{um',k\^) 1 

where the self-energy S is given by (fT2|) . The dependence on system parameters Uq and 7 
can be eliminated if we measure u and T in units of coq, and measure /cy in units of (cjot)^'''- 
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Introducing u = u/uq, T = T/uq, and fcy = k\^/(uojy/'^, we re- write as 

. -■tr^ f dku ^(Urn',k'u) 1 

A;,,) = 2n'XTj2 / " x p _^ ^ . --, (26) 

where, as before, A = l/{2n), and 



f:{oOm,k) = {\/\^\ + ^\-\k\) (27) 

In this section we again consider A not as a given constant but as a free parameter, and 
analyze the limit A ^ 1. The question we address is whether Tc is non-zero already at 
arbitrary small A, or it only emerges when A exceeds a certain threshold. 

Because A is treated as small parameter, Tc is expected to be small, and to get the pairing 
we need to explore logarithmical behavior which comes from frequency scale larger than T^. 
Accordingly, we replace T^^, by (l/27r) f dum' and set ±itT as the lower limits of the 
integration over positive and negative Um', respectively. We then have, instead of fl26l) 

_ - f ^ f ^Ml ^{^m',k\\) 1 

$(a;^, /c||) = ttA / du^> / X - —, (28) 

In general, ^{oJm, k\\) is a function of both arguments, but in proper limits the dependence 
on one of the arguments is stronger than on the other. In the main text we consider the 
limits fcy ^ Cdm and fcy ^ Um- In the first case, the momentum dependence is stronger 
than frequency dependence, and ^{uJrn,k\\) can be approximated by ^{k\\). The fermionic 
self-energy at fcjj ^ Um behaves as 

S(c.„.,^||)^^ (29) 
Substituting this form into (!28|) we obtain 



Integrating explicitly over Um and introducing new variable x = k!^ /T, we obtain Eq 6 in 
the main text. 

$(y)=2A/^^^l^$(x), (31) 

In the opposite limit, fc| ^ w, the momentum dependence of the self-energy is small, 
and S can be approximated by its value at the hot spot, i.e., S ~ \/|^wJ- This limit can 
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only be justified at small or large ratio of Vy/v^ because otherwise substituting this S into 
fl28|) we obtain that typical internal k'^ are of order u . Nevertheless, if we assume that the 
momentum dependence of the fermionic self-energy can be neglected, at least for order-of- 
magnitude estimates, we obtain that ^{um, k\\) can be approximated by and the gap 

equation becomes js 9| 



TT 



(32) 



This is Eq. (9) in the main text. 



2. Solution of the Eq. 

We first replace the soft cutoff imposed by (2VxT + with a hard cutoff, 

^1/^ dx 



= 2A / logx $(x) (33) 

Ji X -\- y 

The integration over x ^ ?/ gives rise to log^ terms in the perturbation theory, which, as we 
know from the analysis in the previous section, do not give rise to the pairing instability. We 
therefore focus on the range x ~ The contribution from this range gives rise to a single 
logarithm (logl/T) if we momentary assume that ^(x) is a constant. Like we did before, 
we treat A as a small parameter and check whether Eq. fl33l) has a solution at a finite T. 

For small A, T is expected to be also small, i.e., the upper limit of the integration in ( 133|) 
is a large number. We first consider y and x ~ ?/ in the range 1 <^ x ^ 1/^, and search 
for the solution of ( l33l) in the form $(x) = exp[— /(p(x))], where p{x) = 27rAlogx. Plugging 
this function back into the equation and introducing a new variable z = x/y, we obtain, 

e-/(2^Aiog,) ^ 2Alogy r ( 1 + ^1 e-«2-Mi°sJ/+i°s-)) (34) 

Ji z + l\ logy J 

y 

We assume and then verify that typical z are of order one, i.e., log 2; <^ logy. Taylor 
expanding in log z / log y we obtain 

""^ dz 1 



1=2 A logy 
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where we introduced 



Solving f p5l) we obtain 



Q{p) 



dm 

dp 



Q{2T[\\ogy) = — arcsin(27rA logy). 



TT 



Integrating over logy, we have 

fip) = 



which is Eq 7 of the main text. 



27r2A 



p arcsm 



in (p) + a/1 — p2 — 1 




(36) 



(37) 



(38) 



FIG. 5: Numerical solution of Q{p) and its comparison with theory. In this figure we have set 
A = 3.7 X 10-^. 



Lastly we consider the boundary condition aX y = = 1/T. For such y, the upper limit 
is 1. Substituting formally the trial solution <l>(x) = exp[— f {p{x))] into the r.h.s. of (133|) we 
obtain 

dz 1 



1 = 2Alog?/o 



(39) 



^0 z + 1 2'9(2'f-'*i°gyo) 

One can easily make sure that to satisfy this equation, Q(27rA logyo) must be larger than if 
we take our result, Eq. (!37|l . and set y = yo there. At vanishingly small T, the presence of 
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the upper limit of the integration over x in (133|) becomes relevant only for y infinitesimally 
close to Ho (for smaller y, the upper limit can be safely set to infinity). This implies that at 
T ^ 0, Q{p) undergoes a finite jump at y = yo, hence at the actual T = 



This last equation is satisfied if 27rAlogl/Tc = 1 Restoring the parameters, we then obtain 
Eq 8 in the main text. For this Tc, the jump in Q at the boundary is between Q = 0.5 (Eq. 
(!37|) at y = yo = 1/Tc) and Q = 0.73, which is the solution of (!39l) at Tc. 

In Fig. [5] we plot Q{p) obtained from the numerical solution of Eq. f l33|) and compare 
it with the (approximate) analytical solution presented in this section. We see that the 
agreement is quite good. 



dQjp) 
dp 



p=27rAlogl/r^ 



= oo. 



(40) 
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